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Abstract 

In this paper, we investigate the stability of a feedback interconnection between a linear system and 
a Duhem hysteresis operator, where the linear system and the Duhem hysteresis operator satisfy either 
the counter-clockwise (CCW) or clockwise (CW) input-output dynamics. More precisely, we present 
sufficient conditions for the stability of the interconnected system that depend on the CW or CCW 
properties of the linear system and the Duhem operator. Based on these results we introduce a control 
design methodology for stabilizing a linear plant with a hysteretic actuator or sensor without requiring 
precise information on the hysteresis operator. 

I. Introduction 

Hysteresis is a common phenomenon that is present in diverse systems, such as piezo-actuator, 
ferromagnetic material and mechanical systems. For describing hysteresis phenomena, several 
hysteresis models have been proposed in the literature, see, for example, flU, [fT8l . [[161 . These 
include backlash model [27) which is used to describe gear trains, Preisach model for modeling 
the ferromagnetic systems and elastic-plastic model which is used to study mechanical friction 
E), |[T8l . From the perspective of input-output behavior, the hysteresis phenomena can exhibit 
counterclockwise (CCW) input-output (I/O) dynamics [[I), clockwise (CW) I/O dynamics [|20l . or 
even more complex I/O map (such as, butterfly map [|3)). For example, backlash model generates 
CCW hysteresis loops, elastic-plastic model generates CW hysteresis loops and Preisach model 
can generate CCW, CW or butterfly hysteresis loops depending on the weight of the hysterons 
which are used in the Preisach model. 



B. Jayawardhana and Ruiyue Ouyang are with the Dept. Discrete Technology and Production Automation, University of 
Groningen, Groningen 9747AG, The Netherlands e-mail: bayujw@ieee.org, r.ouyang@rug.nl 



The CCW and CW I/O dynamics of a system can also be related to certain dissipation 
inequalities flTJ , ll23l . 11261 . Denoting AC as the class of absolutely continuous functions, we 



show in IfTTTl that for a class of Duhem hysteresis operator $ : AC (I 
have that for every u<$, e AC(M + ) there exists a function H 
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for almost every t, where y$ = &(u$,y$, ) and y$ G 1 is the initial condition. The inequality 
© characterizes the CCW I/O property of the operator $. We will discuss this property in detail 
in Section UB Here, we use the symbol O in H Q to indicate the counterclockwise behavior of $. 

As a dual result to IfTTll . in Il24l we give sufficient conditions on the Duhem hysteresis operator 
such that it exhibits CW input-output dynamics. In particular for a class of Duhem operator <3>, 
we construct a function H : IR 2 — » M + which satisfies 
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for almost every t. Correspondingly, the symbol O in H (J indicates the clockwise behavior of 
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Figure 1. Feedback interconnection between a linear plant P and a Duhem operator <I>. 



In this paper, we exploit our knowledge on H and to study the stability of an intercon- 
nected system as shown in Figure [Q where P is a linear system and $ is the hysteresis operator. 
We consider four cases of interconnections where the plant P and the hysteresis operator $ can 
assume either CCW or CW I/O dynamics. These four cases are summarized in Table H 

In Theorem 14. II of this paper, the interconnection case @ in Table U is considered, where both 
the linear system P and the hysteresis operator $ have CCW I/O dynamics. In particular, we 
give sufficient conditions on P which are dependent on the underlying anhysteresis function of 
$ that ensure the stability of the closed-loop system with a positive-feedback interconnection. 



Table I 

Four Possible cases of interconnection 




piczo ) 



This result is motivated by recent results on the positive-feedback interconnection of negative 
imaginary system ll23l and of CCW systems ll26l . A motivating example for the interconnection 
case (a) is the piezo-actuated stages which are commonly used in the high-precision positioning 
mechanisms, see, for example lfl5l . The piezo-actuated stage contains two parts: a piezo-actuator 
and a positioning mechanism, which can be described by 

p mx + bx + kx = I 
V = 

$ : Fpiezo = HV), 

where m is the mass, b is the damping constant, k is the spring constant, c is the proportional 
gain, V is the input voltage of the piezo-actuator, F piczo denotes the force generated by the 
piezo-actuator and x denotes the displacement of the stage. The piezoelectric actuator has been 
shown to have CCW hysteresis loops from the input voltage to the output generated force (see, 
for example [7]). It can be checked that the linear mass-damper-spring system P is also CCW 



(3) 



from F piezo to x or, equivalently, P is a negative-imaginary system 112311 . 

In Theorem I4.3L we consider the interconnection case @ in Table H where the linear system 
P has CW I/O dynamics and the hysteresis operator <fr has CCW I/O dynamics. In this case 
Theorem 14.31 provides sufficient conditions on P which are independent of $ such that the 
closed-loop system with a negative feedback interconnection is stable. An example for this case 
is the active vibration mechanism using piezo-actuator, which has been used for vibration control 
in mechanical structures [fl"3~l . The mechanism can be described by 
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: 



mx, + bx + kx 
V 



F ■ 

A piezo 



piezo ) 
-CX, 



(4) 



As described before, the piezoelectric actuator has CCW I/O dynamics and it can be checked 
that the mass-damper- spring system P is CW from F piczo to x. 

Theorem 15.11 deals with the interconnection case ©, where P has CCW I/O dynamics and 
<fr has CW I/O dynamics. A motivating example for this case is the mechanical systems with 
friction ED . which is given by 

P : mx + kx = -Ff rictioil , 

> (5) 

$ : ^friction = $0), J 

where F ivictioQ is the friction force. As discussed in [1271 . the friction force has CW I/O dynamics 
where the input is the displacement. On the other hand, the mechanical system is CCW from 
the friction force — -Friction to the displacement x. 

As a completion to the Table HI we present the analysis of the interconnection case @ in 
Theorem 15.31 Based on these results, we present in Section |VI] a control design methodology 
for a linear plant with a hysteretic actuator/sensor <fr and we provide two numerical examples 
in Section VII. 

II. PRELIMINARIES 

In this section we give the definitions of the CCW and CW dynamics based on the work 
by Angeli HI and Padthe ll20ll . Figure [2] illustrates the CCW and CW input-output dynamics 
of a (nonlinear) operator G : u \-t G(u) =: y. We denote A(7(R + ,M n ) the space of absolutely 
continuous function / : K + — > M. n . 
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Figure 2. A graphical illustration of counter-clockwise (CCW) and clockwise (CW) I/O dynamics of an operator G : u i — > y. 
(a) CCW I/O dynamics; (6) CW I/O dynamics. 



A. Counterclockwise dynamics 

Definition 2.1: itHl, ||20ll A (nonlinear) map G : AC(M + ,M m ) AC(M + ,M m ) is counter- 
clockwise (CCW) if for every u G AC(M + ,M m ) with the corresponding output map y := Gu, 
the following inequality holds 

liminf / (y(t), u(t))dt > -oo. (6) 

T->oo J Q 

For an operator G, inequality © holds if there exists a function V : IR 2 — > M + such that for 
every input signal u, the inequality 

holds for almost every t where the output signal y := Gu. 

Definition 2.2: A (nonlinear) map G : AC(IR + ,R m ) ->■ AC(IR + ,R m ) is strictly counterclock- 
wise (S-CCW) (see also Q), if for every input u G v4C(IR + , R m ), there exists a constant e > 
such that the inequality 

liminf/" (y(t), u(t)) - e\\y(t) fdt > -oo, (8) 

holds where y := Cm. 

Note that for systems described by the state space representation as follows: 

^ _ x = f(x,u), x(0) = x 1 
2/ = h(x), J 
where x(t) G M n is the state, u(t) G IR m is the input, y(t) G M m is the output and /, h are 
sufficiently smooth functions, the following lemma provides sufficient conditions for £ to be 
CCW (and S-CCW). 

Lemma 2.3: Consider the state space system £ as in ©. If there exists V : R n — > M + and 
£ > 0, such that 



dV(x) r . , /dh{x) r/ . 



dh(x) 



dx 



f(x,u) 
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holds for all x G M n and u G M m , then E is CCW. Moreover if e > 0, it is S-CCW. 



B. Clockwise dynamics 

Dual to the concept of counterclockwise I/O dynamics, the notion of clockwise I/O dynamics 
can be defined as follows. 

Definition 2.4: flU A (nonlinear) map G : AC(R + ,R m ) ->■ AC(R + ,M m ) is clockwise (CW) 
if for every input u E AC(M + , IR' m ) with the corresponding output map y := Gu, the following 
inequality holds: 



For a nonlinear operator G, inequality (flOl ) holds if there exists a function : IR 2 — > IR + such 
that for every input signal u E AC(R + , M. m ), the inequality 



holds for a.e. t where the output signal y := Gu. 

Lemma 2.5: Consider the state space system S as in ©. If there exist a, V : IR m+n — >• M. 
such that V is positive definite and proper, and 




(10) 



dV(y(t),u(t)) 



<(y(t),u(t)) 



(ll) 



dt 



dV(w,x) dV(w,x) 



< (h(x),w) — a(w, x), 



(12) 



dw dx 



f{x,w) 



holds for all x G W 1 , w E R m and q E M m , then E is CW. 

Proof: Define the extended state space system © as follows 



w = q, 



x = f(x,w), 
y = h(x). 

Note that w defines the input in ©. It follows from CL2]) and (Q3]) that 



(13) 



V < (h(x), q) — a(x, w), 



(y,w) - a(x,w) 



which completes our proof by taking w = u. 



III. Duhem Hysteresis operator 
The Duhem operator $ : AC{R + ) xl ->■ AC(R + ), j/$ ) t-)- $(«$, y$ ) =: ?/$ is described 

by 

= fi(y*(t),u<i,(t))u<i> + (t) + / 2 (y*(t),u*(t))u*_(t), y$(0) = y$ , (14) 

where u$+(t) := max{0, «$(*)}, := min{0, u$(i)} and /i : R 2 -)■ R, / 2 : R 2 ->■ R are 

C 1 . We refer to lfl8l . Ifl9l , [|28l for standard properties of the Duhem operator, such as causality, 
monotonicity and rate-independency. 

The existence of solutions to ([141) has been reviewed in [fT8l . In particular, if for every v G R, 
the functions fi and f 2 satisfy 

(7i-72)[A(7i, , y) - fi{l2,v)\ < Ai(f)(7i -7 2 ) 2 , (15) 

(Ti -72) [.Mti ,v) - /2(72, «)] > -A 2 (v)(ti ~72) 2 , 

for all 71, 7 2 G R, where Ai and A 2 are nonnegative, then (fT41 has a unique global solution and 
$ maps AC(R+) xl4 AC(R + ). 

A. Duhem operator with CCW characterization 

To show the CCW properties of the Duhem operator, we review our previous results in [flD . 
In ifm . we define a function H : R 2 — >■ R + for the Duhem operator $ such that © holds 
(under certain conditions on f\ and / 2 ). Before we can define the function Hq for $, we need 
to define three functions which depend on f x and / 2 . 

Firstly, we define a traversing function cu§ which describes the possible trajectory of $ when a 
monotone increasing and a monotone decreasing u<$, is applied to $ from an initial condition. 

For every pair (y$ , it$ ) G R 2 , let y$ , tt$ ) : [tt$ , 00) — > R be the solution of 

z i v ) -y$o = / h{z{a),a) da, Vu G [</*„, 00), 
and let tu$ i2 (-, y$ , w$ ) : (— oo,m$ ] — > R be the solution of 

ziy) -y$ = / f 2 (z(a),a) dcr, G (-oo,u* ]. 



Using the above definitions, for every pair (y$ , m$ ) £ R 2 , the traversing function y$ , w$ ) : 
R — > R is defined by the concatenation of a;$ i2 (-, y$ , w$ ) and w$ ) i(-,2/$ ,«$ ): 

, . ) w#, 2 (v,3/* » w *o) G (-00, u 4o ) 

ws(u,j/* ,tt* ) = < (16) 
I a; 4)1 (u,y$ ,u$ ) Vt;G [tt$ ,oo). 

Again, we remark that the curve u;$(-, y$ , w$ ) is the (unique) hysteresis curve where the curve 

defined in (— oo,m$ ] is obtained by applying a monotone decreasing u§ £ v4C(R + ,R m ) to 

$(m<i>,?/$ ) with m$(0) = m$ and lim^oo w$(t) = —00 and, similarly, the curve defined in 

[m$ , 00) is produced by introducing a monotone increasing n$ £ v4C(R + ,R m ) to 

with «$(0) = m$ and lim^oo = 00. 

The second function we need to define is the anhysteresis function f an , which represents the 
curve where fi(f an (v),v) = f 2 (f an (v),v). 

Another function that is needed for defining H Q is the intersecting function between the 
anhysteresis function f an and the function w$ as defined above. The function f2 : R 2 — > R 
is the CCW intersecting function if w$(fi(7, v), 7, v) = / an (f2(7, vj) for all (7, u) £ R 2 and 
^(7,^) > u whenever 7 > f an (v) and f2(7,t>) < t> otherwise. For simplicity, we assume that 
n is differentiable. In ifTTl Lemma 3.1] sufficient conditions on f x and f 2 which guarantee the 
existence of such are f an be monotone increasing and 

fi(j,v) < ^1^— - e whenever 7 > f an (v) (17) 
dv 

f 2 (rf,v) < ^1^- - e whenever 7 < f an (v) (18) 
dv 

hold with e > 0. 

Theorem 3.1: Consider the Duhem hysteresis operator $ defined in (PT4l) with C 1 functions 
fi, f 2 : R 2 — > R+- Let / an be the corresponding anhysteresis function which is monotone 
increasing and satisfies (fTTT ) and (fTSl . Denote by f2 the corresponding CCW intersecting function. 
Suppose that for all (7, v) in R 2 , ^(7, v) > / 2 (7, f) whenever 7 < / an (i>) and /i(7, v) < f 2 (j, v) 
otherwise. Then $ is CCW with the function H : R 2 — >■ R + be given by 

H h,v) = 'yv- u$(c,7,v) dcr+ w*(<7, 7, v) - / an (cr) da. (19) 

Jo Jo 

Proof: The proof follows from Lemma 3.1 and Theorem 3.3 in [11]. In particular, it is 
shown in 0TJ that 

dgo(y f t W(i)) <fa(*wm). <20) 



where y$ := $(w$,y$ ) and H is non-negative. By integrating (|20l ) from to T we have 



H o(y^(T),u^(T)) - H o (y^(0),u<^(0)) = / y$(r)u$(r)dr. 

Jo 

Since H is nonnegative then 

/ y<j>(T)u<j>(r)dT > -Ha(y$(0),u$(0)) > -oo. 
Jo 



An example of the CCW hysteresis phenomenon is the magnetic hysteresis in ferromagnetic 
material, which has CCW behavior from the input (an applied electrical field) to the output (the 
magnetization). The magnetic hysteresis can be modeled by the Coleman-Hodgdon model Q 
given by 

m(t) = C a \ u$ (t) | [/(u* (*) ) - y* (*)] + u* (t)g(u* (t)), (21) 

where C a is a positive constant, / : R — > R is a monotone increasing C 1 function, such that 
/(0) = and g is locally Lipschitz. The Coleman-Hodgdon model in (1211) can be rewritten into 
the form of ([141) where: 



h(y$, ««) = C Q [/(u$) - y$] + flf(««), / 2 (y*, «*) = -C a [/(«*) - y$] + (22) 

In this case, it has the same structure as in (fT4l) with / on = /. Figure |3] shows the behaviour of 
the Coleman-Hodgdon model using the functions / and g given by 

/(u*) = #0$) = a, (23) 

where b > and a > 0. It can be easily checked that for every u$(t) G R, f\ and / 2 satisfy 
(fT5l) . i.e., for every u$ G AC{R + ) and for every ?/$(0) G R, the solution of ([22]) exists for all 
t G R+. 

Calculating the curve we have 

6a + + (y*(t) - 6u$(t) + ^) e"^-"*) a G [u*(t), oo), 



w*(<t, y*(t),it*(t)) 



6(7 + fef + ~ 6n<1>( ^ + ^0 e c -( CT - n *) tr G (-oo, u$(t)]. 



(24) 



The CCW intersecting function Q(y$(t), «$(£)) is given by 



b — a 

Cq: 

!/*(t)-&u*(t)+^ 

a — b 

y$W-bn*(t)+^ 



V*(t) > fan{u^{t)), 
y$(t) < /an 



(25) 




Since fx and f 2 satisfy the hypotheses in Theorem I3T1 $ is CCW. Denoting u%(t) = £l(y$(t),u<s>(t)), 
we can compute explicitly H Q in (TT9b as follows 



H {y*(t),u*(t)) 



' us,(t)y (t) - ^bu^t) 2 + ^K(t) - u*(t)) 

+£(y*(*) - &u*(*) + fefX 1 - e^C-W-^K*))) y *(*) > f an {u*(t)), 
w*(t)3/*(t) - |6it*(t) 2 + ^K(t) - u*(t)) 

[ - + ^)( e c-(«;(*)-«*(*)) _i) y# (t) < f an ( U0 (t)). 



(26) 



The graphical interpretation of #0 is shown in Figure 0] where the value of at a given 
time t is given by the area in grey. 

Proposition 3.2: Consider the Duhem operator $ satisfying the hypotheses in Theorem 13.11 
Suppose that / an (0) = 0. Then the function Hq(',v) (where H is as in CH])) is radially 
unbounded for every v. 

Proof: Let us consider v > 0. To show the properness of H (-,v), let us first consider the 
case where 7 > f an (v). In this case, we rewrite the function H , as follows 



7 - /anO)d(X + 



n(-y,v) 



Due to the property of the CCW intersecting function £7, 7 > f an (v) implies that ^(7, v) > v. 

Hence the last term on the RHS of the above equation is non-negative, i.e., J^ 1 '^ o;$(cr, 7, v) — f an (cr)da > 





H (y q> (t),uM)) _ 

L ^4^/* 




yr^ fan 


' J U<f>(t) 

fi(y$(t),M$(t)) 



Figure 4. Graphical interpretation of the function H,^(y<s>(t),u<s,(t)) of the Coleman-Hodgdon model using / and g as in {23} 
with b = 5 x 1(T 3 , C« = lx W 2 , a = 2.5 x 1(T 3 and y* = 0. 



0. Then, 

PV PV 

H (-y,v)> / 7 - fan(<r)do- > / 7 - / an (f )drr = (7 - c)v , (27) 
Jo Jo 

where c := f an {v). Equation (ITTT) indicates that for every t> > 0, Hq{^,v) — >■ 00 as 7 — >■ 00. 

To evaluate the other limit when 7 — >■ — 00, let us consider the case when 7 < 0. Note that 
in this case 7 < f an (v) due to the monotonicity assumption on f an and / an (0) = 0. Rewriting 
H , we have 

pQ(-f,v) pv 

H o{l,v)= 7 - f an (a)da + 7 - u)$(a, 7, v)der 

</o Jn{i,v) 

n(f,v) r0 

7 - ian0)d<7 = / /an(0") - 7d<7. 

Jn(-y,v) 

The last inequality is obtained due to the property of the CCW intersecting function tt, where 
Q(j,v) < v whenever 7 < f an {v). Since ujq, is monotone and non-decreasing (due to the 
positivity of fx and f 2 ) and using the fact that f an is monotone increasing and f an (0) = 0, it 
can be checked that 7 < implies that ^(7, v) < 0. 

Now let us fix 7 such that > 7 > 7- Using the fact that a>$(cr, 7, v) > u$(cr,j,v) for 
all a < v and using monotonicity of f an , it follows that > > ^(7, v) where the constant 




n := 0(7, v). Thus 



H (l,v) > / /an(o-) - 7da 

7n(7,u) 



= / /an(c) - 7d<7 + / / arl ( CT )- 7d (7 

Jsi(f,v) Jn 

> / /an(o-) - 7^0- > / /an(n) - 7d<T 



= (7-/ an (fi))a 

The last equality shows that as 7 — > —00, iJ — > 00 since < 0. Therefore, we can conclude 
that for the case v > 0, the function Hq(-,v) is radially unbounded. 

Using similar arguments we can get the same conclusion for the case when v < 0. ■ 

B. Duhem operator with CW characterization 

The CW property of the Duhem operator has been discussed in our previous results in [|24l . 
where we also constructed a function H : IR 2 — > IR + for the Duhem operator such that © 
holds. Following a similar procedure as before, the construction of the function H requires 
three functions: the traversing function the anhysteresis function f an and the intersecting 
function A. The definitions of the functions uq, and f an are the same as those given in Section 
IIII-AL However the CW intersecting function A has a different definition than that of the function 

a 

The function A : R 2 — > R is a CW intersecting function if o;$(A(7, v ), 7, v) = / an (A(7, v)) 
for all (7, v) G R 2 and A(7,t>) < v whenever 7 > f an (v) and A(7,t>) > v otherwise. Here we 
assume A is differentiable. In ll24l Lemma 1] sufficient conditions on fi and / 2 which ensure 
that such A exists are f an be monotone increasing and 

/l ( 7)?; ) > ^M + e whenever 7 > f an (v) (28) 

f 2 (j, v) > d ^7^ + e whenever 7 < f an (v) (29) 
df 

hold with e > 0. 

We recall our main results in [24] in the following theorem, which gives the sufficient 
conditions on $ such that it is CW. 



Theorem 3.3: 11241 Theorem 1] Consider the Duhem hysteresis operator $ defined in ([141) 
with C 1 functions fx, f 2 : — > M + . Let f an be the corresponding anhysteresis function which 
satisfies (1281) and (1291) . Denote by A the corresponding CW intersecting function. Suppose that 

for all (j,v) in M 2 , fi(j,v) > f 2 (j,v) whenever 7 < f an (v) and fi(j,v) < f 2 {l,v) otherwise. 
Let the anhysteresis function f an satisfies f an (0) = 0. Then $ is CW with the storage function 



H : M 2 — »■ R + be given by 



fan{v)d<J 



ui$(a, 7, f )dcr, 



(30) 



The proof is similar to that of Theorem 13.11 where we use also the result in [24, Theorem 1] 
which shows that H satisfies ©. 

An example of the CW hysteresis phenomenon is the friction-induced hysteresis in mechanical 
system, which has CW behavior from the input (i.e., the relative displacement) to the output 
(i.e., the friction force). One of the standard model to describe friction-induced hysteresis is the 
Dahl model ||6l, ETI . which is given by 



sgn(ii$(t)) 



sgn(u$(i)) u$(t) 



(3D 



sgn \1 

where y§ denotes the friction force, u$ denotes the relative displacement, F c > denotes the 
Coulomb friction force, p > denotes the rest stiffness and r > is a parameter that determines 
the shape of the hysteresis loops. 

The Dahl model can be described by the Duhem hysteresis operator (fT4l with 



1 _ m 



sgn 1 - — , f 2 (y®,u< s> ) = p 



y<s> 

1+ k 



In Figure [51 we illustrate the behavior of the Dahl model where F c = 0.75, p = 1.5 and r = 1. 

It is immediate to check that fx and f 2 satisfy the hypotheses in (fT5l) . which means that for 
all m$ G AC(M + ) and y<s>(0) G R the solution of (I3TT) exists for all t G IR+. The anhysteresis 
function of the Dahl model is f a n(u$(t)) = 0. 
Calculating the curve we have 

F c + (y*(t) - F c ) e £(«*(*>-"> a G [«*(*), 00), 
-F c + (yz(t) + F c ) et^W) a G (-00, u*(t)]. 
The CW intersecting function A(y$(t),u<s>(t)) is given by 

y*{t)>o, 
y*{t) < 0, 



o;$(cr, 



(33) 



A(y*(t), «*(*)) 



u$(t) + —In- ... ,. 



u*(t) - fin 



(34) 



V*(t)-F c 



t 10 "- 6 u{t) 

(a) (b) 
Figure 5. The input-output dynamics of the Dahl model with F c — 0.75, p = 1.5 and r = 1. 




Figure 6. Graphical interpretation of the function H(j(y&(i),u$(t)) of the Dahl model using fi and fa as in J32b with a = 1, 
F c = 0.5 and y$ = 0. 



Denoting = A(y$(t), tt$(t)), we can compute explicitly the function H (J in (1301) as 

follows 

' -F e (u*(f) - «&(«)) + f (y*(t) + F c )(l - e £W*>-«*W)) y s (t) > 0, 



i?o(y* (*),«*(*)) 



The graphical interpretation of H is shown in Figure |6j where the value of H CJ at a given 
time t is given by the area in grey. 

Proposition 3.4: Consider a Duhem operator $ satisfying the hypotheses in Theorem 13.31 
Suppose that f an is monotone increasing and / an (0) = 0. Then the function H(j(-,v) (where 



is as in (1301) ) is radially unbounded for every v. 

Proof: To show the properness of H (-,v) for any given v, we first consider the case 
7 > fan{v). Since the Duhem operator $ satisfies the hypotheses in Theorem 13 .31 the function 
H (J is nonnegative. Thus, using (1301) and since A(f an (v),v) = v we have 

#o(7, «) > #0(7, - H (f an (v),v) 

rMl> V ) / , A(7,«) pA(f an (v),v) 

= / /an(^)dcr- / w$((j,7,t;)da+ / w*(«r, 7, u)d(7 

J Mfan( v ) > v ) J v J v 

/•A(7,«) 

= / / an (cr) -a;$(cr,7,f)dcr. 

By the definition of the CW intersecting function A, 7 > f an {v) implies that A (7, u) < t>. Using 
the mono tonicity of w$, w$(cr, 7, i>) < 7 for all cr < t>, and thus, it follows from the above 
inequality that 

/■A(7,u) 

H {-y,v) > / /an(o-) -a; # ((7,7,u)d(7 

/•A(7,w) 

> / /an(o-) - 7dcr, 

Now let us fix 7 s.t. f an (v) < 7 < 7. Since u;$(er, 7, w) < u>$(<7, 7, v) for all a < v and using 
the monotonicity of f an , we have that A(7,i>) < A, where A = A(j,v). Therefore, 

/•A(-f,v) pA 
Ho(l,v) > / fani?) ~ 7do" > / / an ((7)-7d(7 

V J V 



A 



> / /an(u)-7d<r = (7-c)(t;-A) >0 

where c := f an (v) and A = A(7,f) for any given v. Hence, it implies that for every v, 
Hoili f ) — c» as 7 — t- 00. 

We can apply similar arguments to show that for every v, #0(7, u) — > 00 as 7 — >■ —00 by 
evaluating the case when 7 < f an (v). ■ 

IV. Linear system with CCW Duhem hysteresis 

In this section we analyze the stability of a feedback interconnection of a linear system 
and a CCW Duhem hysteresis operator. The stability of the closed-loop system is analyzed by 
exploiting the CCW or CW properties of each subsystem. 



(35) 



Theorem 4.1: Consider a positive feedback interconnection of a minimal single-input single- 
output linear system and a Duhem operator $ as shown in Figure Q] satisfying the hypotheses 
in Theorem 13.11 as follows 

x = Ax + Bu, 

P : 

y = Cx, 

* : = /i(y* (*)>«*(*) )"*+(*) + /ifa* (*),«*(*))"*-(*)> 

where A G M nxn , B G M nxl and C G M lxn . Let £ := (CB)- 1 where we assume that CB > 0. 
Suppose that there exist £ > and Q = <5 T > such that 

- (A T Q + QA) + eA T C T CA < 0, (36) 

QB + A T C T = 0, (37) 

Q - iC T C > 0, (38) 

hold and the anhysteresis function f an satisfies {f an {v) — £,v)v < for all v G M. (i.e. f an belongs 
to the sector [0, £]). Then for every initial condition (x(0),y$(0)), the state trajectory of the 
closed-loop system (1351) is bounded and converges to the largest invariant set in {(x, y$)\CAx + 
CBy„ = 0}. 

Proof: Using V(x) = \x T Qx and (1361) and (|37|). it can be checked that 

V = l - x T {A T Q + QA)a; + x T QBu 
< -ex T A T C T CAx - x T A T C T u 

= (u T B T C T + x T A T C T )u - e(CAx + CBuf(CAx + CBu) 

= (y,u) -zy 2 - 

It follows from Lemma 12.31 that the linear system is S-CCW. 

By the assumptions of the theorem, the Duhem operator $ is also CCW with the storage 
function H ■ IR 2 -> R + as given in (fT9l) . 

Now let H c [(x, y$) = V(a;) + H (y$, Cx) — Cxy<$, be the Lyapunov function of the intercon- 
nected system (|35l . We show first that if c / is lower bounded. Substituting the representation of 



V and H a , we have 

fCx /•Q(y^,,Cx) 



H c i = -x T Qx + zCx — / w$(o", Ca;) da + / u;$(cr, y$, Cx)da 

2 Jo Jo 

- / / on ((r)d(7 - Cxy$ 

Jo 

= -X T Qx- / /an(0 r )d(T+ / W$ (d, ?/$ , Cx) - f an (a)da 

^ JO J Cx 

j /-Cx /-Co; i>£l(y$,Cx) 

>-x T Qx- / (/an(cr) - £cr)dcr - / £ada + / u;$(cr,y $ ,Ca;) - / an (cr)dcr 

^ JO JO JCx 

i rSi{yi>,Cx) 
> -x T (Q-£C T C)x + / u*(a,y*,Cx) - f an (a)da. (39) 

^ JCx 

where the last inequality is due to the sector condition on f an . In the following, we will prove 
that the last term on the RHS of ([39]) is lower bounded. Notice that since f x > 0, f 2 > 0, (fTTT ) 
and (fT8l) imply that ^"'^ > e for some e > 0. Hence f an is strictly increasing and invertible. 

Consider the case when y§ > f an (Cx) which implies also that VL(y$,,Cx) > fan(y<s>) > Cx 
by the definition of fi. Using the monotonicity of uq, we have 

rtt(y$,Cx) ^{y^,Cx) rfmiy®) 

j u$(a, y$, Cx) - f an (v)da > y$ - f an (a)da > y$ - f an (a)da. 

JCx JCx JCx 

Define V(y*, Cx) := f*£ {m) - f an (a)da and let c := ^" /a 2 " (Ca) + f an (Cx). It follows that 
fan'M > /"» and f an (a) < c for all a e [Cx, /-»]. Therefore 

V(y$, Cx) > yg,- fan{°)de > y$ - c da 

JCx JCx 
= l(y*-fan(Cx))(f- n 1 (c)-Cx)>0. 

Thus, ^ _ f an (a)da is lower bounded by V(y<s>, Cx) which is positive definite (it is 

equal to zero only if y$ = f an (Cx)) and V(y$, Cx) — > oo as yq, — > oo. 

When y$ < f an {Cx), we can obtain the same result where J^ y9 ' Cx > y^ _ / an (cr)dcr is lower 
bounded by I^(|/$,Cx) which is positive definite and V{y$,Cx) — > oo as y$ — > — oo. 

Therefore, using (l39l , we have 

#d > \x T {Q - iC T C)x + V(y 9 , Cx), 
which is radially unbounded. 



Now computing the time derivative of H d , we obtain 

H d = V + H - Cxy® - Cxy$ < -ey 2 . 

This inequality together with the radially unboundedness of H c i imply that the trajectory (x, is 
bounded. Using the Lasalle's invariance principle, we conclude that the trajectory (x, of (1351) 
converges to the largest invariant set contained in M := {(x, G R n x R|CAc + CBy§ = 0}. 

■ 

We illustrate Theorem 14.11 in the following simple example. 
Example 4.2: Consider 

P : x = —x + u, y = x, 

* : m = fi(y<s>(t),u$(t))u<s, + (t) + f 2 (y$(t),u<s,(t))u^(t), 
u = y$, u$ = y, 

where x(t) G R and the functions fx, f% satisfy the hypotheses in Theorem 13.11 Using Q — 1, it 
can be checked that d36l) — d38l) hold. Using H d as in the proof of Theorem 14.11 let us define 
H c i(x,y$) = \x 2 + H (y$,y) — yy$ and a routine computation shows that 

H c i < yy$> + ${y)y - yy$ - yy^ - y 2 
= -(-x + y®) 2 . 

Note that Q = 1, C = 1, so that (l50l) holds for ^ < 1. This means that the result in Theorem 
14. II holds if the anhysteresis function f an satisfies (f an (v) — £v)v < 0, for all v G R and £ < 1. 
In other words, / a „ should belong to the sector [0, £] for the stability of the closed-loop system. 

A 

The result in Theorem 14.11 deals with a positive feedback interconnection of a linear system 
and a Duhem hysteresis operator. This is motivated by the study of an interconnection between 
counterclockwise systems as studied in 0] for the general case and in [|23l for the linear case. 
In the following result, we consider the other case where a negative feedback is used instead. 

Theorem 4.3: Consider a negative feedback interconnection of a minimal single-input single- 
output linear system and a Duhem operator $ as shown in Figure \T\ satisfying the hypotheses 



in Theorem 13.11 as follows 

x = Ax + Bu, 

P : 

y = Cx + Du, 



(40) 



* : = /i(y*(*),«*(*))w*+(*) + f2(y<5>(t),u$(t))u0-(t), 

u = -y®, uq> = y, 

where A G M nxn , B G R nxl , (7 G M lxn and Del. Assume that there exist P = P T > 0, L 
and 5 > such that the following linear matrix inequalities (LMI) 



P [ »xi] - [c*-] > 
J(^[rT] + [A %]P)+6L T L<0, 



(41) 
(42) 



hold. Then for every initial condition (x(0), j/$(0)), the state trajectory of the closed-loop system 
(|4Q|) is bounded and converges to the largest invariant set in {(x,y$)\L = 0}. 

Proof: By the assumptions of the theorem, the Duhem operator $ is CCW with the function 
H : R 2 — > JR + as given in (fl9i 

Define the extended state space of the linear system in (l40l) by 

w = q, 

Pext '■ X = Ax + BW, 



(43) 



where w = u. 

Using V = \[w x T ] T P 



w 

X 



y = Cx + Dw, 



, a routine computation shows that 



1, 



w x T ] ( 


B T 




o nxn 




w 




gnxn J^T 


P +P 


B A 


) 


X 



V 



Using d4T]) and (|42]> . 



+ [ w x T ]P 



1 

gnxl 



(44) 



This inequality ((44)) with q = u (by the relation in (l43l ) implies that the linear system defined 
in W is CW. 



Now take H d (x,y^) = H a (y$,Cx — Dy$) + V(x,y§) as the Lyapunov function of the 
interconnected system (|40t . where H d is radially unbounded by the non-negativity of H and 
the properness of V. It is straightforward to see that 

H d = H + V, 

< (v,u) + - 5 

= -5 









2 




-y* 






L 








X 













2 




-2/* 






L 








X 







(45) 

where the last equation is due to the interconnection conditions u = — y® and y = u§. It follows 
from (|45T ) and from the radial unboundedness (or properness) of H c i, the signals x and y$ are 
bounded. 

Based on the Lasalle's invariance principle ifTTl . the semiflow (x, y§) of (|4Q|) converges to the 
largest invariant set contained in M := {(x,y<s>) GR"x K|L [~f*] = 0}. ■ 
To illustrate Theorem 14.31 let us consider the following simple example. 
Example 4.4: Let 

x — ~~ \~ 
y = —2x + u, 

u = = i/, 

where x(t) E M and the functions fx, f% satisfy the hypotheses in Theorem 13.11 By using 
P = [ } 2 ~q 2 ] , it can be checked that (|4T|) — ((42)) hold. Following the same construction as in 
the proof of Theorem @7; 
computation shows that 



P 



the proof of Theorem 14.31 we define H d (x, y$) = \x T Px + H Q (—y$, —2x + y$) and a routine 



#d < -2(-3z + y$) 2 + yy* - $(y)y 



-2f-3x 



2/$ J 



Thus, we can conclude that (x, y$) converges to the invariant set where x = 



A 



(46) 



V. Linear system with CW Duhem hysteresis 

Dual to the result that we present in the previous section, the feedback interconnection of a 
linear system and a CW Duhem hysteresis is considered in this section. 

Theorem 5.1: Consider a negative feedback interconnection of a minimal single-input single- 
output linear system and a Duhem operator $ as shown in Figure \T\ satisfying the hypotheses 
in Theorem 13.31 as follows 

x = Ax + Bu, 

P : 

y = Cx, 

* : m = /i(j/*(f),u* (*))«*+(*) + /2(j/« (*),«*(*))«»-(*), 
u = -y®, u$ = y, 

where A e R nxn , B e M nxl and C e M lxn . Let e := {CB)- 1 where we assume CB > and 
assume that there exist Q = Q T > such that 

- {A T Q + QA) + £A T C T C71 < 0, (47) 

Qfi + A T C T = 0, (48) 

hold. Then for every initial condition (x(0), y$(0)), the state trajectory of the closed-loop system 
(|46|) is bounded and converges to the largest invariant set in {(x, y$)\CAx — CBy§ = 0}. 

Proof: Let V(a;) = \x T Qx, and using (1471)— (1481). it can be checked that 

V = ^x T (A T Q + QA)x + x T QBu 
< (y,u) - ey 2 . 

It follows from Lemma 12.31 that the linear system is S-CCW. 

By the assumptions of the theorem, the Duhem operator $ is CW with the function H CJ : 
R 2 -»■ K+ as given in (|30b . 

Now let H c i(x,y<$>) = V(x) + H Cj (yq>,Cx) as the Lyapunov function of the system (l46l) . 
According to Proposition I3.4L H (y<f>, Cx) is radially unbounded for every x, which implies that 
Hd(x,y$) is radially unbounded. 

Computing the time derivative of H c i, we obtain 

H cl = V + H < -ey 2 . 



This inequality together with the radially unboundedness of H d imply that the trajectory (x, y<$,) is 
bounded. Using the Lasalle's invariance principle, we conclude that the trajectory (x, y$) of (|46l) 
converges to the largest invariant set contained in M := {(x, G M n x R|CAc — CBy§ = 0}. 

■ 

To illustrate Theorem 15.11 we could use the same linear system as given in the Example 14.21 
Example 5.2: 

P : x = —x + u, y = x, 

® m = /i(2/*(f), «*(*))«*+(*) + / 2 (2/*(t), «*(*))«*_(*), 

u = -y$, uq, = y, 

where x(t) £ R and the functions fi, f 2 satisfy the hypotheses in Theorem 13.31 Using Q = 1, 
it can be checked that (|47T) and (1481) hold. Define H c i(x,y<s>) = \x 2 + H (y^,,y), sl routine 
computation shows that 

H d < yy$ - yy$ - y 2 
= -(-x + y$) 2 , 

which implies that (x, y§) converges to the invariant set where x = y$. 

A 

Theorem 5.3: Consider a positive feedback interconnection of a minimal single-input single- 
output linear system and a Duhem operator $ as shown in Figure \T\ satisfying the hypotheses 
in Theorem 13.31 as follows 

x = Ax + Bu, 

P : 

y = Cx + Du, 
® m = fi(y<s>(t),u$(t))u$, + (t) + f 2 (y<f>(t),u$(t))u<s,-(t), 
u = y$, u$ = y, 

where A e W lXn , B e R nxl , C e M lxn and Del. Assume that there exist P, L, 5 and rj > 
such that P = P T > rj [ qt c ] > and the following linear matrix inequalities (LMI) 

P^xi] = [ qt ] j (50) 
U P [b°T} + [ ^ B a t]P)+SL t L<0, (51) 



(49) 



hold. Assume further that /i(7, v) < \ and ^(7, v) < | for all (7,1*) G M 2 . Then for every 
initial condition (x(0),y§(0)), the state trajectory of the closed-loop system (|49l ) is bounded and 
converges to the largest invariant set in {(x, yq>)\L [ y g] = 0}. 

Proof: Define an extended system P ext as in (l43l) and let V(w,x) = \ [w x T ] P Using 
(|50l). (|5TT) and (03]), we have 



V<(y,u)-5 



L 















(52) 



Equation (|52|) indicates that the linear system is CW. 

Next we take H d (x,y<s>) = H (yq,,y) + V(yq >) x) — yy<$> as the Lyapunov function of the 
interconnected system. We will show first that H d is lower bounded. Using the definition of V 
as above and H (J as in (|3Q|) . we have 

<"A(?/ S ,u 4 >) 



H, 



1 r 











T 

if or 


P 










X 





/ an (cr) - w$((T, 2/$, M.j,)dcr+ / w$(cr, «$)da-n$y$ 



D 2 DC 



Since uq> = y (by the interconnection), w| = [ w x T ][ i 
there exists rj, e > such that P — rj [J^ D qTq\ > zl- Then 



j [ 1 ] . By the assumption on P, 



H, 



cl 







D 2 DC 




w 


T 

W X 


(p-v 


C T D C T C 


) 








X 



V 2 



A(y$,«*) 



f an (a) - u;$(er, w$)dcr 



U<J, 



+ 



e 

> - 
~ 2 



o;$(cr, y$, it$)dcr — u<$>y<$> 
2 

+ 



J /»-^ lW) «,)d ff + y ^w*(«7,y*,u«) - y$ + -w$J do- 



(53) 

It can be checked that the second term of (1531) is nonnegative. Indeed, it follows from the property 
of the CW intersecting function A that if A(y$, > we have that / an (c) > cu$(cr, 7/$, u$) 
for all a" G [0, A(y$,u$)] and if A(y$,M$) < then f an (a) < w$(cr, y§, u<&) for all a G 
[A(j/*,m*),0]. 

To check whether the last term of (|53l is lower bounded, we use the definition of u<$, given 
in the Section [Tll-Al Consider the case m$ > 0. Using the definition of tu<$> in (fT6l) . the last term 



of ( 1531 can be written by 

J (w*(<r,y#,u*) - + do- 

y$+ 72(^0(5, s)ds ) dcr + 



f/<|, 



- - / 2 (o;$(s, tt$), s)dsd(x + ~u| > 0, 

where the last inequality is due to fact that / 2 (7, v) < \ for all (7, v) G M 2 . In a similar way, 
we can obtain the non-negativity of J^* (cj$(cr, w$) — ?/$ + |w$)dcr for the case < 0. 
Therefore, <153T > implies that H d is lower bounded and radially unbounded. 
It can be computed that 



Hd = V + H - yy$ - yy$ < -5 



L 










X 





(54) 



P 



Hence, by the radially unboundedness of H d , (|54l) implies that (x, y$) is bounded. Using the 
Lasalle's invariance principle, we can conclude that the trajectory (x, y<&) of (l49l) converges to 
the largest invariant set contained in M := {(x, y<s>) 6R"x W\L [ y * ] = 0}. ■ 

To illustrate Theorem 15.31 let us consider the Example 14.41 where we replace the negative 
feedback interconnection by a positive one. 

Example 5.4: 

x = —3x + y$, 
y = -2x + y$, 

y$ = fi{y${t),u$(t))u$ + (t) + /2(y*(*),u* (*))«*-(*) 

u = uq, = y, 

and the functions fi, f 2 satisfy the hypotheses in Theorem 13.31 By using 
P = [} 2 ~ 6 2 ] , the conditions in (EH and (ED hold with L = [1 -3] and 5 = 2. Also P = P T > 
\ [-2 ~4 2 ]' i- e - V=i- Usin S H d (x,y$) = \ [j» x T ] P [ y £] + H (y$, -2x + y 9 ) - yy^, routine 
computation shows that 

Hd < -2(-3x + ^) 2 . 

Hence, if fi(j,v) < \ and f-z{"j,v) < \ for all (7, v) G M 2 , then (x,y$) converges to the 
invariant set where x = —^y<s> following Theorem 15.31 

A 



where x(t) G 



G 



G 



C 



(a) 



(b) 



Figure 7. Feedback interconnection of a linear plant G, controller C and hysteresis operator <&. (a) An interconnection example 
where the plant G is driven by a hysteretic actuator 4>; (b) An interconnection example where the dynamics of G is measured 
by a hysteretic sensor <I>. 



VI. Controller design 

The stability analysis given in the previous sections can be used to design a controller for a 
linear plant with hysteretic sensor/actuator. Consider the closed-loop system as shown in Figure 
[71 where G and C are the linear plant and controller, respectively, and they are given by 

G I x G = A G x G + B G u Gl c I x c = A c x c + B c u G , 
( Vg = C G x G + D G u G , [ y c = C G x G + D c u G . 

Thus depending on the location of the hysteretic element, the cascaded linear systems can be 
compactly written into 

x = Ax + Bu 

(56) 

y = Cx + Du, 

where x = [%c \ an d for the case of hysteretic actuator as shown in Figure Ha), A = [ B ^ c a c ] » 
B = [ Be Dal' C = [ DcCg g °]> D = D G D G , or for the case of hysteretic sensor as shown in 
Figure Hb), A = [ A G B( f c c ] , B = [ Bg b ° c c ],C=[c a d g c c ], D = D G D C . The controller design 
can then be carried out as follows. 

• Control design algorithm for the case of CCW $: 

1) Determine the anhysteresis function f an of the Duhem operator $ and possibly, the 
desired L. 

2) Find C such that either (|36l> - (|38ll or d4TT>-(l42l holds. 

3) If (f36])-(f38l) is solvable, then C stabilizes the closed-loop system with a negative 
feedback interconnection; otherwise 

4) If (|4TI)-(14"2~1) is solvable, then C stabilizes the closed-loop system with a positive 
feedback interconnection. 



• Control design algorithm for the case of CW $: 

1) Determine the functions fi and f 2 of the Duhem operator $ and possibly, the desired 
L. 

2) Find C such that either (@7])-(|48]> or (|50]>-(|5D holds. 

3) If (I47l)-(l48l) is solvable, then C stabilizes the closed-loop system with a negative 
feedback interconnection; otherwise 

4) If (l50l)- ((5n) is solvable, then C stabilizes the closed-loop system with a positive 
feedback interconnection. 

Putting (|56l) into the setting of our main results in Theorem 14.31 14.1[ 15.11 and 15.31 the invariant 



set is contained in M := {(xq, xc, y$)\N 
design parameter for determining C. 



XC 



0} where the matrix iV can also become a 



VII. Numerical examples 



x 



m 



b 



Figure 8. Mass-damper-spring system connected with a hysteretic actuator 



As an example, we consider a mass-damper-spring system with a hysteretic actuator denoted 
by as shown in Figure [8l where m is the mass, b is the damping constant, k is the spring 
constant and x denotes the displacement of the mass. Let m = 1, b = 2 and k = 1, then the 
mass -damper- spring system is given by 





f 1 \ 






X = 






: 




l-l -2 1 







V = 1 )x + u. (57) 



A. CCW hysteretic actuator 

Let us first consider the case when the hysteretic actuator has CCW I/O dynamics, such as 
piezo-actuators [15]. Assume that the actuator is represented by the Duhem operator (fl~4)) where 

^(7, v ) = - 7 + 0.475v + 0.3, / 2 (7, v) = 7 - 0.475v + 0.3, V(7, v) E R 2 . (58) 

It can be verified that f an (v) = 0.475-u and the functions fx and / 2 satisfy the hypotheses given 
in Theorem 13.11 

With A c = [i 2 i 4 ], 5 C = [0], C c = [-1.5-2] and D c = 1, conditions (gB-dH are 
- 1 1 -1.5 -2 1 



solvable with P 



and L = [0 1/4 0]. Hence the controller 



1 7.74 5.51 -8.74 -15.86 
5.51 7.4 -5.51 -14.36 
-1.5 -8.74 -5.51 10.24 17.86 
_ -2 -15.86 -14.36 17.86 38.36 _ 

C can stabilize the closed-loop system with negative feedback interconnection. In this case, 
iV = [0 1/40 0]. According to Theorem 14.31 the velocity of the mass -damper- spring system 
converges to zero and the position of the mass-damper-spring system converges to a constant. The 
closed-loop system is simulated in Matlab/Simulink with the initial condition x(0) = [—10 5] T 
and the results are shown in Figure 0a). 

On the other hand, since we have f an ( v ) = 0.475t>, then by taking A c = [_° 2 i 4 ], B c = 



oi C c = [i l] and D c = 0, it can be checked that (I36h-(l37l) holds with f = 0.5 and Q = 

1 -6 -2 "| 

In this case iV = [10-2-31]. Moreover, f an belongs to the sector [0,0.5]. 



1 4-1-4 

-6-17 3 
^-2 -4 3 7 _ 

Similar to the previous case, it follows from Theorem 14. II that the velocity of the mass-damper- 



spring system converges to zero and the position of the mass-damper-spring system converges 
to a constant. The simulation results is shown in Figure Hfb). 

B. CW hysteretic actuator 

For the case of a CW hysteretic actuator, see for example the magnetorheological (MR) damper 
used in the structure control [25], the mass -damper- spring system is given by (|57[) . Assume that 
the actuator is represented by the Duhem operator (Q~4) where 

/i(7,«)= 0.25(1 -7), f 2 (7,v) = 0.25(1 + 7), V( 7) w)6l 2 (59) 

The anhysteresis function for this Duhem operator is f an = 0. It can be shown that fx < 0.25 
and f 2 < 0.25 for all v 6 R. In addition fx and / 2 satisfy the hypotheses in Theorem 13.31 hence 
the Duhem operator with (l59l is CW. 



6 



10 



(a) 



(b) 



Figure 9. Simulation results of the numerical example with CCW hysteretic actuator, (a) The negative feedback interconnection 
case with the initial condition :r(0) = [—10 5] T ; (b) The positive feedback interconnection case with the initial condition 
x{0) = [-10 10] T . 



With A c 
with P = 



1 ], B c 



= [Q], C c = [i l] and D c = 0, the conditions (|47T> -(|48 l) are solvable 
Hence the controller C can stabilize the closed-loop system with 



-2 -4 J 
'5 1-5-2' 

1 3-1-3 
-5-16 3 
_ -2 -3 3 6 _ 

negative feedback interconnection. In this case, N = [1 — 2 — 3 1]. According to Theorem 



15.11 the velocity of the mass -damper- spring system converges to zero and the position of the 
mass-damper-spring system converges to a constant. The simulation results are shown in Figure 
[TOf a) with the initial condition x(0) = [10 5] T . 

Since we have f x < 0.25 and f 2 < 0.25 for all v G R, by taking A c = [_° 2 J 3 ], B c = [§], 
C c = [-3-i] and D c = 2, it can be checked that ((50|) -(|5l ]) holds with 6 = 1, ri = 0.5, 

r 2 2 -3-11 

2 3n.8fi 15.83 -32.8R -2fi.9 

L = 01/4000 and P 



2, it can be checked that (I50l)-(I5TI) holds with 5 

2 2 -3 -1 

2 30.86 15.83 -32.86 -26.9 
15.83 38.26 -15.83 -51.4 
-3 -32.86 -15.83 35.86 27.9 
_-l -26.9 -51.4 27.9 74.54 _ 

velocity of the mass-damper- spring system converges to zero and the position of the mass 



. It follows from Theorem 15.31 that the 

to zero and the position of the mass- 
damper-spring system converges to a constant. The simulation results is shown in Figure [TOl b). 



VIII. Conclusions 

It has been shown in this paper that the stability analysis of a linear system with hysteresis 
nonlinearity is accommodated by exploiting the I/O property of the corresponding hysteresis 
operator. Furthermore, the stability analysis enables a straightforward control design methodology 
for a plant with hysteresis nonlinearity without having to know precisely the parameters of the 
hysteresis operator. It offers a different paradigm in the design of controller for such systems 




20 40 60 80 100 20 40 60 80 100 

time s time s 

(a) (b) 

Figure 10. Simulation results of the numerical example with CW hysteretic actuator, (a) The negative feedback interconnection 
case with the initial condition x(0) = [10 5] T ; (b) The positive feedback interconnection case with the initial condition 
x(0) = [10 - 5] T . 



where we do not need to define an inverse hysteresis operator which is commonly used in 
practice. The dissipativity approach which is used in this paper can be extended directly to 
nonlinear plants with hysteresis nonlinearity. One possible class of nonlinear plants which can 
be treated with our approach is the CCW systems as studied by Angeli [1J and by van der Schaft 
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